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Abstract 

We compute the tensor CMB anisotropy power spectrum for singular and 
non-singular instantons describing the beginning of an open universe accord- 
ing to the Euclidean no boundary proposal. Singular instantons occur gener- 
ically, whereas non-singular instantons require more contrived scalar field po- 
tentials. For the latter, we consider potentials in which a sharp feature, either 
negative or positive, is added to a gently sloping potential. In the first case one 
finds a nearly divergent contribution to the low multipole CMB anisotropy, 
in conflict with the COBE observations. In the second case the divergence is 
weaker, but matching the low multipoles forces the added feature to be large 
and narrow. For singular instantons, there is a better match to the obser- 
vations, without any such contrivance. The distinction between singular and 
nonsingular instantons disappears in the limit as the universe becomes flat, 
but is still observable for densities as high as 0.7 of the critical density. 

I. INTRODUCTION 

In the most common approach to inflationary theory one postulates a scalar field with a 
gently sloping potential, and assumes that for some reason the field was initially displaced 
from the potential minimum. If the initial displacement is large, the field approaches a slowly 
rolling state in which the universe inflates. This state is an attractor, and in it the system 
loses memory of the initial conditions. This scenario, which is certainly the simplest version 
of inflationary theory, predicts that the universe should be flat to high accuracy today. It also 
predicts that the initial state of the universe should be totally inaccessible to observations 
today, since the scales most relevant to defining the initial state were stretched by inflation to 
scales currently exponentially larger than the Hubble radius. If future measurements confirm 
the universe is very nearly flat, then, assuming inflation is the explanation, discussions of 
what came before inflation although interesting will remain strictly academic. 
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Current CMB observations are consistent with a flat universe, for example the recent 
Boomerang measurement ^ yields 0.65 < Qtot < 1-45 at 95 per cent confldence. This lends 
support to the hope that the simplest version of inflation, described above, might be correct. 
However, signiflcant space curvature is not yet excluded by the observations. This paper 
is devoted to examining the observational consequences of inflationary scenarios in which 
signiflcant space curvature would exist today, and in which the initial conditions for the open 
universe are actually visible in the microwave sky. In an open universe, the curvature scale 
of the universe on the surface of last scattering subtends an angular scale of approximately 



i/n^ radians, about 25 degrees for f2o = 0.3. If we live in such a universe, cosmic microwave 

sky observations can probe the initial conditions for the inflating universe. 

Theories of open inflation were initially constructed from scalar fleld potentials with false 
vacua, using instantons known as Coleman-De Luccia instantons [Q. These can be inter- 
preted as describing tunnelling from a prior false vacuum inflationary state [§,0, although 
the relevant instantons only exist for rather special potentials. More recently, however, it 
was realised that open inflation can occur far more generically through a class of singular, 
but flnite action, instantons which exist for essentially all gently sloping inflationary po- 
tentials. The regular instantons do have the virtue that the prediction of Q is unique in a 
given theory. For the Hawking- Turok instantons, the most probable universe a priori is one 
with a very low value of Q, but there are solutions for essentially all values of Q up to unity. 
In the absence of a better understanding of how the actual value is determined, which may 
involve some sortof anthropic considerations, we shall here simply treat the value of f2 as a 
parameter to be adjusted to flt the universe we see. The pattern of density perturbations is 
then, for given fl and given scalar potential, uniquely predicted. 

In this paper we exhibit an interesting observable difference between non-singular and 
singular instantons. We discuss a generic problem faced by non-singular instantons and 
show how it is alleviated in singular instantons. 

II. GRAVITATIONAL INSTANTONS AND OPEN INFLATION 

Instantons are saddle point solutions of the Euclidean path integral, and open inflationary 
instantons may be naturally interpreted within the framework of Euclidean quantum gravity 
and the no boundary proposal 0. The instantons provide a saddle point, which one can 
expand around to compute the Euclidean path integral. Correlators of interest are uniquely 
deflned in the Euclidean region, and then analytically continued into the Lorentzian universe. 
We have recently carried this program through to leading (quadratic) order for scalar and 
tensor perturbations (Related calculations were performed in a different approach 

in refs. [jTT|,|12[.) The well known problems of the non-positivity of the Euclidean Einstein 
action and the non-renormalisability of quantum gravity do not enter in these low order 
calculations. 

Until recently the class of known cosmological instantons was quite limited. Coleman 
and De Luccia discovered the flrst examples when generalising the problem of the decay of 
a false vacuum in scalar fleld theory to include gravity [||. In the limit of weak gravity the 
decay is well understood and occurs via bubble nucleation. In a localised region of space the 
scalar fleld quantum tunnels through the barrier stabilising the false vacuum. The bubble 
so formed expands at the speed of light and inside it the scalar fleld rolls down to the true 
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FIG. 1. Inflationary potentials of the types being considered here. The left is a smooth function, 
like m?(j)^ , A(/>^ or e*^*^. In this theory only Hawking- Turok singular instantons exist. On the right 
are two potentials allowing Coleman-De Luccia instantons. The solid line shows a potential with 
a sharp false minimum added, the dashed line one with a sharp maximum. 

vacuum. In the presence of gravity, instantons only exist for scalar field potentials with a 
sufficiently sharp false vacuum (as shown in Figure 1). The reason is that the gravitational 
instanton has finite size, ~ Mpi/ \/V where Mpi is the Planck mass and V the potential 
energy density. For an instanton to exist, in which the scalar field is not constant, the 
scale of variation of the field must be smaller than the instanton size. But this scale of 
variation is determined by the second derivative of the potential in the region of the barrier, 
|V^0(^| = M^. The condition for existence of Coleman-De Luccia instantons is therefore that 

Coleman-De Luccia instantons may be used to describe the nucleation of bubbles in a 
false vacuum region of de Sitter space 0,0] • The interior of such bubbles then form infinite 
open universes and with modest fine tuning of the distance A0 over which the field rolls 
during inflation, one can adjust the value of VLtot to an interesting value less than unity 
today. But in order for the Coleman-De Luccia instanton to exist, the condition mentioned 
above must be satisfied. Assume for example that the potential is approximated by |m^(/)^ 
in the neighbourhood of the true vacuum, where m « Mpi. For N efolds of inflationary 
expansion, one requires to roll for 2\/NM^i where Mpi is the reduced Planck mass. The 
false vacuum has to be at least this far from the true vacuum. But existence of the Coleman- 
De Luccia solution requires >> 4iVm^, and for reasonable iV > 40 (for acceptable Vl 
today), the scale M must be at least an order of magnitude larger than m. As we show 
later, yet another tuning is required in order to suppress the large angle CMB anistropies. 

In open inflation, it was assumed that the scalar field became stuck in the false vacuum, 
leading to large amounts of inflation, in the course of which the universe approached perfect 
de Sitter space. Bubbles would nucleate in this de Sitter space, as the field tunnelled 
through the barrier between the false and true vacuum states. Each bubble contains an 
infinite, inflating open universe 

Coleman-De Luccia instantons have a different interpretation in the no boundary pro- 
posal. There they are viewed as classical solutions describing the rounding off of the 
Lorentzian universe on a compact Euclidean region. The Euclidean path integral uniquely 
specifies the spectrum of fluctuations inside the bubble without the need for additional as- 
sumptions regarding the pre-bubble era. Indeed all calculations that have been performed 
to date have in effect used the instanton background to define the pre-bubble era. This 
makes the calculations identical to those performed in the no boundary interpretation. To 
that extent one can say that all predictions of open inflation are really predictions of the no 
boundary proposal. The pre-bubble inflating universe appears to be a redundant theoretical 
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construction. 

The Coleman-De Luccia instantons are interesting because they provide a calculable 
scenario for open inflation. As mentioned the potentials needed to obtain such instantons 
are necessarily contrived (Figure 1). However Hawking and one of us recently showed that 
a class of singular but finite action Euclidean instantons exists for almost every gently 
sloping inflationary potential We have computed the spectrum of fluctuations about such 
singular instantons and found that in spite of the singularity the correlators are uniquely 
defined, just as in the Coleman-De Luccia case 0,^. Interesting differences arise because the 
singularity imposes Dirichlet boundary conditions on the perturbation modes. We pointed 
out that the observational effect of this difference is likely to be most pronounced in the 
tensor spectrum and this is what we discuss here. We shall show that the part of parameter 
space for Coleman-De Luccia theories in which the bubble size is much smaller than the De 
Sitter radius, so that the tunnelling is very similar to that in flat spacetime, is ruled out. 

III. FLUCTUATIONS ABOUT SINGULAR AND NON-SINGULAR 

INSTANTONS 

As we shall see, Coleman-De Luccia instantons with potentials as shown in Figure 1 
generically produce a very large amplitude of long wavelength tensor modes. Let us make 
clear at the outset however that this constraint cannot be used to rule out all such models. 
The tensor perturbations are governed by the height of the inflationary potential, and can be 
adjusted independently of the scalar perturbations by flattening the potential. For example 
potentials of the form used in hybrid inflation, with a flat plateau followed by a sharp drop 
produce acceptable scalar perturbations but almost no tensor component. Another way of 
suppressing the low multipoles is to further tune the potential so that it is steep around 
the values of where the bubble nucleates 0. In the light of our discussion above, this 
constitutes a third fine tuning, needed to make such models work. 

In this paper we compare the predictions of singular and Coleman-De Luccia instantons 
with potentials of the form shown in Figure 1. In the singular case we assume a simple 
monomial potential like (jP' or 0^, and in the Coleman-De Luccia case we superpose a sharp 
negative false vacuum. We show that unless Vttot is rather close to unity today the Coleman- 
De Luccia examples are generally ruled out because they predict unacceptable large angle 
anisotropics in the microwave sky. In contrast the singular instantons which occur generically 
in gently sloping inflationary potentials appear more compatible with the observations. 

In previous work we have derived the two-point correlators of the scalar and tensor |^ 
metric perturbations in open inflationary universes associated with both classes of Euclidean 
cosmological instantons. All perturbations are determined from correlators of the gauge- 
invariant Newtonian potential and the transverse traceless tensor perturbation tjj, which 
may be computed directly from the path integral. 

To first order in ^ the Euclidean correlators are specified by a Gaussian integral P,^. For 
both regular Coleman-De Luccia instantons and singular instantons the result is unique. In 
the latter case the singularity enforces Dirichlet boundary conditions. The Euclidean two- 
point correlators are analytically continued into the Lorentzian region where they describe 
the quantum mechanical vacuum fluctuations of the various fields in the state described by 
the no boundary proposal initial conditions. 



4 



In the present work we shall calculate the temperature fluctuations on the microwave 
sky from the Lorentzian two-point correlators. The key observable difference between the 
two types of instantons occurs for wavelengths of order the curvature scale. Since the long- 
wavelength continuum in the scalar power spectrum vanishes linearly with wavenumber p, 
the differences are small there. Likewise the bound state of the scalar perturbation potential, 
producing long range correlations beyond the curvature scale in the open universe, is known 
to have a very minor effect on the CMB anisotropy |T^. However the spectrum of primordial 
gravity waves has for the regular instantons a near divergence at small p and therefore 
provides a better opportunity for a distinguishing test. 

The result for the symmetrised two-point correlator tensor metric perturbation about 
either Hawking-Turok or Coleman-De Luccia instantons is 

{{U,{x)A'A^')}) = fcothpvr^,(r)^_,(rO +r,^^^M^) ^^jf^ (1) 

Jo p \ smhpvr J a[Tja{T') 

where k = SnG, and length units are chosen so that the spatial curvature scale is unity. In 
this formula r is the conformal time as defined in and x the comoving radial coordinate. 
The bitensor W^^jfy (x) is the sum of normalised rank- two tensor harmonics with eigenvalue 
= + 3) of the Laplacian on |Q. The eigenmodes 5'p(r) are solutions of the 

Lorentzian tensor perturbation equation 

/ a" \ 

+ - - 1 gp{r) = p^gpir) (2) 



y dr'^ a 

normalised to obey gpir) —>■ e"*^"^ as r ^ — oo. 

First note the potential divergence in the integrand due to the 1 / (psinhpvr) in the 
first term. This divergence, as we shall now argue, is cancelled by the second term. The 
second term involves Vp which is a reflection amplitude computed in the Euclidean region. In 
conformal coordinates the metric takes the form b'^{X){dX'^ + dQ^). For singular instantons 
we have < X < oo where the singularity is located at X = 0. For regular instantons 
we have — cxd < X < cxd. In both cases the perturbations obey a Schrodinger-like equation 
with potential U{X) = — 1. This potential diverges to +oo at X = in the singular 
case, but is finite everywhere in the regular case. In fact in the latter case it is close to a 
refiectionless potential —2 sech^X. The quantity Vp is in both cases the reflection amplitude 
for waves incident from X = -|-oo. For singular instantons it is by unitarity a phase but for 
non-singular instantons it is a complex number of modulus less than unity, and it is small 
at high p. Both reflection amplitudes tend to minus one as p — because long-wavelength 
modes are completely reflected, hereby yielding an infrared finite correlator. However, since 
the non-singular Coleman-De Luccia instantons are much closer to the perfect 5*^ non- 
reflecting solution, we expect the reflection amplitude to tend to —1 at much lower p than 
in the singular Hawking-Turok case. 

The region of low p in the the tensor spectrum is what is known in the literature as 
the bubble wall fluctuation spectrum [IU|. When the de Sitter symmetry is only weakly 



broken, with a scalar field present, there is a near divergence in this spectrum. These long- 
wavelength tensor perturbations give a substantial contribution to the CMB anisotropies. 
From the discussion above, we expect a larger contribution to the large angle microwave 
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anisotropies for regular instantons. In other words, the mild breaking of de Sitter invariance 
in non-singular models allows for large long-wavelength fluctuations about the background 
solution. On the contrary, in singular models the deviation from an 0(5) instanton is drastic, 
and the singularity keeps such fluctuations small. 

IV. NON-SINGULAR "THIN- WALL" INSTANTONS 

The scalar field equation in the Euclidean region reads 

{b'<P,^),^ = b'V^, (3) 

where a is the proper radial distance {da = hdX). Following p we consider the case where 
the potential is given by superimposing a sharp negative "bump" of amplitude —AV centred 
about (f)f onto a smooth monotonically increasing function of 0. On a non-singular instanton, 
the scalar field rolls in the upside down potential from 0o; gaining kinetic energy until it hits 
the "bump" and rapidly decelerates to an almost standstill near (j)f. Effectively all of the 
kinetic energy of the field is converted to potential energy and any damping is negligible. 
The field then remains approximately constant as the scale factor h turns round and vanishes 
as (um — c")- 00 is fixed by the form of the potential and the requirement of regularity. This 
generally implies that the scalar field must have reached the "bump" well before the scale 
factor turns round. We therefore take h ^ a m. equation (||), and approximating V^^ as 
y,(t>o = (0o), we have 

^ crV<^o- (4) 
We have = at the regular pole, so we may solve to find 

~ 00 + ^ V4,y. (5) 

If the field approaches the "bump" at cxfe, then its kinetic energy just before hitting the 
"bump" is 

l€ ^ ^ \y<^o - M (6) 

and we may equate this to AV. 

As we shall discuss shortly, it is useful to rewrite the Schrodinger equation in the Eu- 
clidean region in a form where it involves the potential U = 10'^ where prime denotes 
derivative with respect to conformal coordinate X. The strength of the potential is then 

[ ^<p'^dX = I -h<p\da ^ ^Yi^_ (7) 
J 2^ i 2 ^''^ 2 64 ^ ^ 

If we take the smooth part of the potential to be of the form A0"', we may introduce 
the quantities = ^ and = f ^ (0o)- N is the slow roll approximation to the number 
of inflationary efoldings in the open universe. H is the slow roll Hubble parameter, with 
6 ~ 4t sin Ha. Then we can write 
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9n (Hab) , . 

^ - 256N ■ 

We will see below that for regular instantons the quantity C provides an infrared cutoff for 
the amplitude of the bubble wall fluctuations. From the condition that the scalar field must 
have reached the bump well before the scale factor turns round, using 6 ~ sin Ha we 
see that Hab must certainly be less than |. As a concrete example, if we take n = 2 and 
N = 50, this yields C < 0.01. Generically, in the regime where the bubble radius is much 
smaller than the radius of the de Sitter space, C will be very much smaller than this, since 
the formula involves the fourth power of the size of the bubble, at- 



V. NON-SINGULAR "THICK- WALL" INSTANTONS 

We now consider instantons associated with potentials with a sharp positive feature 
as shown by the dashed curve in Figure 1. In this case, the scalar field motion is confined 
within the region of the feature over the instanton, and does not probe the smooth part of the 
potential at all. Unlike the thin-wall case discussed above, the scalar field varies significantly 
over the whole instanton, and not just over a localised region of it. The starting value of the 
scalar field is tuned so that the field reaches the peak of the feature at approximately the 
same time as the scale factor rolls over. If indeed the potential is exactly symmetrical about 
the peak over the region probed by the instanton, these two events occur at exactly the 
same moment. The change in sign of the slope of the potential may then be able to balance 
the antidamping, bringing the scalar field to a halt as the scale factor again tends to zero, 
giving us a non-singular solution. For this to occur, the feature must be sufficiently sharp. 
This can be achieved using different iable functions with large curvature at the peak. We can 
model this by introducing a kink. We model the potential in the vicinity of the feature at 
say 0=1, as V — 10 — 0*1, with V and constant and positive. Then we approximate b as 
1/ H sin Ha, with H^ = ^V, assuming that V dominates over gradient energy in the field. 
In this approximation reaches 0,, at a = 7t/2H, and 0o- is odd about a = tt/2H. So in 
order to calculate C = J ^hcp^da, we need only work out (p^^ up to cr = 7r/2if and multiply 
by two. From the scalar field equation we have 



V, 



sin^ Hada 

(9) 



sin^ Ha Jo 

Va> (cos^ Ha -3 cos Ha + 2 



SHsin' Ha 



and so 



C = 2x 






2 9if2 Jo 
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~ 4^ 





COS Ha — 3 cos Ha + 2 da 



(10) 

We can also integrate Eq. (|) to find that A0 = 0, - 0(0) = 1/2(1 + 2\n2)V^/{KV) ^ 
1.19V (I,/ (kV). Inserting into ([T0|) we see that C can be expressed two ways, either as 
C ^ A(f)V^/V = AV/V, or as C ^ O.8A0VM|,,, where M|,, = k'^. We have checked that 
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the above expressions match the numerically calculated values quite closely up to C ~ 1. 
In order to get a value of C close to unity, one requires a large feature in the potential - i.e. 
a large change in V to occur over a range of which is at least of order unity in (reduced) 
Planck units. The calculations shown below exclude small values of C, corresponding in the 
thick wall case to small positive features on the potential. 



VI. EUCLIDEAN REFLECTION AMPLITUDES AND MODES 

The primordial gravity wave spectrum is given by equation ([^). In terms of the proper 
distance a we used in the previous section, we shall fix the integration constant involved in 
defining the conformal coordinate X by setting X = -^^^ non-singular instantons 

we follow references and and define at to be that value of a for which b is maximum. 
For singular instantons it is taken instead to be the value of a at the singularity. 

For singular instantons, the singularity acts to impose Dirichlet boundary conditions 
in the Euclidean region. The only allowed mode function at fixed p is given by ipp 
apC^P'^ + a_pe~*^"^ as X — > cxd. In the non-singular case, we have left and right moving 
modes. The left mover is g^f^ {X) e'^P^ , as X ^ -oo and g^^^ (X) Cpe'P^ + dpC'^P^ 
as X ^ +00. These mode functions satisfy the differential equation 

^ +U{X)]gp{X)=p'gp{X). (11) 



V dX^ 

which has a trivial bound state solution b{X) with = —1. This corresponds to a constant 
shift in the metric perturbation which is pure gauge. It is very convenient to project this out, 
since its presence means that there is an extra phase shift of tt produced by the potential 
even at very low p. The projection is simple. Rather than g one considers g = b j which 
is clearly zero for the bound state ||rT|. This variable also satisfies a Schrodinger equation 



+U{X)]-gp{X)=p%{X) (12) 



dX 



where U is the positive-definite quantity mentioned in the previous section. We define 
(jp (X) in an identical fashion. From the constancy of the Wronskian and using b ~ e"'"^' at 
the regular poles, one finds Vp = j^fp. For singular instantons the reflection amplitude Vp 
is given by the phase and in the non-singular case it equals j-. It is straightforward to 
calculate Vp numerically for any background instanton of interest. 

For the non-singular instantons considered in the previous sections tJ {X) is sharply 
peaked around a value of X, Xb say. We can then make a very good analytic approximation 
to as follows. We replace U{X) by the delta function potential C6{X — Xb) of equiv- 
alent strength, with C as defined in Eq. (|^ or Eq. ((0)) as required. We can then solve 
analytically for g^^^^ {X) and find 



dp 1 + S 



(13) 
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We approximate Xi, as follows. In the "thin-wall" case Xt, corresponds to cr = cri,. Then 
with b ^ sin Ha, 

Xb ^ r -^^da ^ - In tan ^ (14) 

where we have used Hat ~ f • the "thick- wall" case U (X) is simply peaked around at 
and so X^ ^ 0. 



VII. TENSOR CMB ANISOTROPY IN OPEN INFLATION 



The Euclidean no boundary path integral allows us to compute correlation functions 
of any observable. If these correlations are well approximated by a classical statistical 
distribution, as macroscopic observables such as the microwave anisotropies are, we can 
regard the predictions as being characterised by the classical distribution. Our observed 
universe is one member of this classical ensemble. To compare different theories with regard 
to an observation carried out in our universe, we compute how likely the given observation 
is according to each theory. 

We consider the microwave background anisotropy generated by primordial fluctuations, 
expanded in spherical harmonics 6T/T = YjCLimXim {,0,4>) ■ The o/m's obey a^^ = ai-m and 
we have 2/ + 1 real observable quantities for each /. Rotational invariance implies that 
the 2Z + 1 quantities are independently distributed with zero mean and common variance 
Cj^. Neglecting higher-order effects, their probability distributions are Gaussian. For a 
given / we average over the squares of the 2Z -|- 1 observable quantities in our universe to 
determine 0°^^. Then for a given theory of this type, Cf^^/Cf^ is x^-distributed over the 
ensemble of universes with 2/ + 1 degrees of freedom. Cf^^ itself is gamma-distributed 



with probabihty density function / (c°bs; + i) /cf^^ ^ + |) El 



If 0^"°' 



the median value of C, 



th 
I ) 



is greater than 

then the fraction of universes with Ci less likely than Cf^^ 



IS 



given by 21 (/ + |, (/ + |) Cf^yCj"^) /T (/ + |). Similarly, if Cf^s is less than the median 



value of C;*^, 



then the fraction of universes with Ci less likely than C"^'^ is given by 2 



2T (/ + i,(/ + i)cr^7c*'^) /r(/ + i 



We need to obtain the C* s for the different theories we are interested in. Using the 



usual Sachs- Wolfe formula this is given in terms of our symmetrised tensor correlator 
(S) as 



Cf' = t^^ dr rfr'fcothpvr $^(r)$^p(r'^ 



1 



sinh pn 



(15) 



The primordial tensor power spectrum at the end of inflation defines inital conditions for 

which 



the Sachs- Wolfe integral. To compute the multipole moments we use CMBFAST |16 



evolves the mode functions from the surface of last scattering at riss up to the present time 
Tnow, given the initial power spectrum. Modifications were required to improve the resolu- 
tion at low wavenumbers, necessary for the accurate evaluation of the the low / multipoles. 
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We then combine the tensor component in the correct ratio with the standard scale 
invariant scalar spectrum of perturbations in order to obtain the total Cf^ to compare with 
experiment. To extract the primordial tensor power spectrum from equation (0), we first 
construct approximate solutions for the eigenmodes (r) = $p (r)a(r). In the inflationary 
phase of the open universe the mode functions closely follow perfect de Sitter evolution in 
which they tend to a constant after the physical wavelength has been stretched outside the 
Hubble radius. Hence, to determine the amplitude and phase of this constant we approx- 
imate a (t) as sinh Ht until fl is close to one, and introduce the associated conformal 
coordinate 

r) = — — ; (it = Intanh — , (16) 

' Jt sinh Ht 2 ' ^ ^ 

7] ^ —oo being the start of inflation, and ?7 — as the universe continues formally to inflate 
without end. r — 77 is a finite constant during inflation whilst this approximation for a{t) is 
a good one. The approximate Lorentzian tensor perturbation equation is then 



+ ^:^ fp(v)=p'fl:{v) (17) 



y dr]"^ cosh rj 
and the solution in which fp{f]) —>■ e"*^'' as 17 ^ —00 is 

^ tp — I 

At a given value of t then, with corresponding r and r], we have gp (r) ^ fp{^)- So 

dividing by a and taking the late-time limit we see that 



-ipTQ 



%{to)^-H, -. (19) 

ip — 1 

Here tq is the conformal time as defined in at the end of inflation. This can be calculated 
numerically and is 0(1) for singular instantons and 0(0.01) for "thin-wall" non-singular 
instantons. From equation (^ the primordial tensor power spectrum Pt{p) at the end of 
inflation is 

2k^- fcothp7r$,(ro)<l>_,(ro) + 1^^- M^lMo)^ . (20) 
p \ 1 + tp smh pir J 

For singular instantons Vp = dp/d_p is a phase factor and can be written as e^*^^. So the 
tensor power spectrum P^{p) for singular instantons is 

in this approximation. For a given potential one evaluates 6p numerically and obtains an 
empirical fit. In the long- wavelength limit dp —>■ — 7r/2 so the power spectrum is infrared 
finite. Actually, it turns out that the CMB power spectrum predicted by singular instantons 
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differs only a little from the one with a perfect reflecting potential in which the ratio Op/ a_p 
is replaced by —1 for all p. 

For non-singular instantons, we have fp = Cp/dp in (^). Using our approximations for 
this in the previous section, we obtain 



P-(p) = (tanh^ + ^ ( 1 - ^"^^^(^^ + -o) + ff^^^(^^ + -o) , , . (22) 

^ p{l+p^)\ 2 sinhpvr ^ 1 + ^ II V ; 



Equations (^if ) and ( ^21) define the initial conditions for the numerical computation of 
the Sachs-Wolfe integral for the different models. A Taylor expansion of the second term 
around p = shows that in the "thin-wall" speculated earlier for typical values 

of C the regime Cp/dp — — 1 sets in at much lower p than in the singular Hawking-Turok 



case. One can see from (|T^) that this leads to a larger contribution to the large angle 
microwave anisotropics for regular "thin-wall" instantons. In the next section we discuss to 
what extent this characteristic feature of false vacua models allows one to observationally 
distinguish them from singular open inflation models. 



VIII. NUMERICAL RESULTS 

In Figure we compare the CMB anisotropy power spectra for singular and nonsingular 
instantons, in open universes with Qt^t = 0.3 and 0.7 respectively. For the nonsingular 
instantons there is a large contribution from the tensor component, shown by the dashed 
line. The amplitude of the large angle contribution is governed by the parameter C discussed 
above. We argued above that on general grounds C has to be smaller than 0.01 for 'thin- 
wall' instantons. For 'thick-wall' instantons C can be larger if the feature on the potential is 
large. For singular instantons there is no such parameter to vary. We have chosen C = 0.025 
for the nonsingular instantons which is certainly conservative for the thin wall case. The 
divergence at low / would be even more pronounced in the allowed regime. 

These calculations show that even if the curvature of the universe today is quite modest, 
one nevertheless sets strong constraints on the form of the inflationary potential and on the 
nature of the primordial instanton. As emphasised above, C decreases as the fourth power 
of the size of the bubble wall thickness - if the latter is much smaller than the Hubble radius 
of the de Sitter space, C is much smaller than 0.01. Since the amplitude of the correlator 
diverges as / dpp~^, the amplitude of the quadrupole diverges roughly as C~^. 

In the panels shown the result for the gravity wave spectrum explained above has been 
combined with the usual scalar spectrum of perturbations appropriate for an open universe 
0. The ratio of tensor to scalar contributons is a function of cosmological and model 
parameters. However, for medium multipoles, (/ ~ 30), the ratio approaches its well known 
flat space value |T^. This value then fixes the relative normalisation of scalar versus tensor 



anisotropy for all multipoles. For a A0" inflaton potential the flat space ratio i?fl = 0.05n. 
Therefore, the higher the value of n is, the more important the contribution from the bubble 
wall tensor fluctuations. In the plots shown we have taken the ratio to be that for n = 2. 
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FIG. 2. Cosmic microwave sky predictions of different cosmological instantons. The upper two 

panels show predictions for an open universe, with fltot = 0.3 and no cosmological constant, for 
Hawking-Turok (panel A) and Coleman-De Luccia instantons (panel B). The former is for an 
^m^0^-potential, the latter for a model where a false vacuum has been added (see text). The lower 
panels compare the Hawking-Turok (C) and Coleman-De Luccia (D) theories for an ^matter = 0-3, 
^Lambda = 0-4, ^tot = 0-7 cosmology. The difference at low / is still marked. These results are for 
a cold dark matter dominated universe with ^cdm = 0.25, baryon density = 0.05, and Hubble 
constant h = 0.65. 
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This yields a quadrupole ratio R2 ~ 0.13 in a singular model and R2 ~ 0.57 in a regular 
model. Higher values for n would allow us to exclude the nonsingular models more strongly. 

For both singular and non-singular instantons we see the rise in Cf^ at low /, character- 
istic of an open universe. We compare these different models to the COBE DMR data as 
follows. First of all have to set the overall normalization of each model. We do this with 
the RADPACK software 
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Using the DMR data alone, we find the normalization 
which maximizes the likelihood for each model. We then compare likelihoods amongst the 
different models. The relative likelihoods follows: 

Singular : Singular with A : Non-singular : Non-singular with A : Flat spectrum 
28 : 76 : 1 : 22 : 97 



where the flat spectrum is one with constant Z (/ + l)C/'s, shown for comparison. In 
this, Bayesian, approach, the nonsingular instantons for an open universe with Qfot = 0.3 
are strongly disfavoured. 

Having done the likelihood analysis above, we now take an hypothesis-testing approach, 
using the probabilities given in Section VII above. The strongest constraint on the mod- 
els comes from the quadrupole, and we focus on that here. Unfortunately the true sky 
quadrupole is not yet known, and the literature contains various estimates of it. We have 
therefore assumed a range of values taken from various references. It is to be hoped that 
the MAP experiment will accurately determine the actual value. 

The model dependence of the low I Ci's suggests that in order to quantify the difference 
between the models we should compare the correctly-normalized predicted quadrupole mo- 
ments with the measured quadrupole moment. In the table below we show the percentage of 
universes in the ensemble associated with a given theory with a measured quadrupole more 
extreme than that seen. We compare the results from singular and non-singular instantons 
with the best fit fiat spectrum for comparison. We have done this for a selection of groups' 
estimates for the observed quadrupole |20|-^. We have converted all measured values to 
the dimensionless quantity / (/ + 1) Q/27r, dividing by (2.73 K)^ where necessary, to match 
the output of CMBFAST. The result from [^] is effectively a direct measurement of the 
quadrupole, albeit with a systematic error due to the galactic cut. The other results are 
harder to interpret, having been obtained using maximum likelihood techniques with highly 
non-Gaussian likelihood functions for the quadrupole ||21|. This means that the quoted val- 
ues below should have large skewed error bars. We also show what a measurement of a 
larger quadrupole (that from the best fit fiat spectrum) would tell us for illustration. 



Measured value 

3C2/71 


Singular 
^]A = 
1.7 X 10-1° 


Singular 
= 0.4 
1.2 X 10-1° 


Non- singular 

nA = o 

4.0 X 10-1° 


Non-singular 

= 0.4 
2.3 X 10-1° 


Flat spectrum 
1.0 X 10-1° |21| 


0.11 X 10-1° 


20| 


0.56% 


1.2% 


0.071% 


0.27% 


1.9% 


0.20 X 10-1° 


21| 


2.3% 


4.9% 


0.31% 


1.1% 


7.0% 


0.37 X 10-1° ^ 


9.0% 


18% 


1.3% 


4.6% 


25% 


1.0 X 10-1° |21 


61% 


96% 


13% 


37% 


83% 



We note that in general the probability is is several times larger for the singular case 
as compared to the the non-singular case. Both models are easier to rule out at a given 
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confidence level than the flat spectrum. Note that even this model is ruled out at the 98% 
level if the result of [20| is taken at face value! It should be remembered that for our non- 
singular "thin-wall" model, we have assumed a value of C even larger than the extreme 
best-case. In situations where the theoretical quadrupole is much larger than the measured 
one, the probability scales as {C™^^/C^^ ~ C^/^ for the non-singular models. The "thin- 
wall" nonsingular instantons appear to be strongly ruled out by the observed quadrupole, 
even if VLtot is as large as 0.7. 



IX. CONCLUSION 

We have computed the tensor CMB anisotropy power spectrum for a class of singular and 
non-singular instantons. We showed that this provides a way to observationally distinguish 
different versions of open inflation. The "thin-wall" false vacuum models generate larger 
fluctuations on large angular scales, distinguishing them from singular models. Using the 
COBE data, we have found that this characteristic feature strongly disfavours "thin-wall" 
Coleman-De Luccia instantons relative to the singular Hawking- Turok models. Non-singular 
"thick-wall" Coleman-De Luccia instantons are still viable, but only if the false vacuum 
feature in the scalar potential is large. In this case the predictions depend strongly on the 
detailed parameters describing the feature and the models are hence somewhat unattractive. 
These calculations have therefore enabled us to further constrain the form of the inflaton 
potential in open inflation. 
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